Periodicity-dependent stiffness of periodic hydrophilic-hydrophobic hetero-polymers 
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From extensive Monte Carlo simulations of a Larson model of perfectly periodic heteropolymers 
(PHP) in water a striking stiffening is observed as the period of the alternating hydrophobic and 
hydrophilic blocks is shortened. At short period and low temperature needle-like conformations 
are the stable conformation. As temperature is increased thermal fluctuations induce kinks and 
bends. At large periods compact oligomeric globules are observed. From the generalized Larson 
prescription, originally developed for modelling surfactant molecules in aqueous solutions, we find 
that the shorter is the period the more stretched is the PHP. This novel effect is expected to stimulate 
polymer synthesis and trigger research on the rheology of aqueous periodic heteropolymer solutions. 
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Almost all the important "molecules of hfe", e.g., 
DNA, RNA and proteins, are hetero-polymers [Q . There- 
fore, in order to gain insight into the in- vivo "structure" 
and "function" of these macromolecules, in recent years 
physicists and chemists have been studying the in-vitro 
structure and dynamics of simpler hetero-polymers con- 
sisting of only two different types of monomers. The se- 
quence distribution is totally random in what are known 
as random heteropolymers (RHP) B. The RHP are of 
special interest to theorists also because of their close re- 
lation with the random energy model 1 1 and spin glasses 
M; these similarities and the unusual properties of the 
RHP are consequence of the combination of quenched dis- 
order and a special type of frustration arising from the 
competing interactions in the RHP MM- Very recently, 
random heteropolymers with correlated sequence distri- 
bution has also been considered theoretically . On the 
other hand, perfectly periodic heteropolymers (PHP) 
have begun to receive attention only very recently [g|. 
Orlandini and Garel |H| carried out what may be loosely 
called the first in-vacuo g] Monte Carlo (MC) simula- 
tions of PHP. The aim of this paper is to report the re- 
sults of in-vitro MC simulations of a very simple model of 
PHP in water to demonstrate a novel dependence of the 
stiffness of the PHP on the periodicity of the hydrophilic 
(or, hydrophobic) segments. 

We follow the recent reformulation [pO| of the Lar- 
son model [|ll|,Q of surfactants in water ||l^ to model 
the PHP in water. In the spirit of lattice gas models, 
the system is modelled as a simple cubic lattice of size 
Lx X Ly X Lz- Each of the molecules of water can oc- 
cupy a single lattice site. A surfactant occupies several 
lattice sites successive pairs of which are connected by a 
nearest-neighbour bond of fixed length. We shall refer to 
each site on the surfactants as a monomer. The primary 
structure of each PHP can be described by the symbol 

IpOpIpOp IpOp where / and O refer to the hydrophilic 

and hydrophobic monomers and the basic building block 
IpOp, each of length Lp, is repeated n times such that 
2Lpn — La is the total length of the PHP. No monomer 



is allowed to occupy a site which is already occupied by 
a water molecule. Besides, no two monomers of the PHP 
are allowed to occupy the same site simultaneously. 

If the chain consisted of only hydrophilic monomers it 
would behave exactly as a self-avoing walk in vacuo be- 
cause of the complete identity between the hydrophilic 
monomers and the molecules of water. On the other 
hand, if it consisted of only hydrophobic monomers it 
would collapse forming a compact globule. What makes 
the model PHP so interesting is the competition be- 
tween these two conformations arising from the compet- 
ing hydrophilic-hydrophobic effects. 

For the convenience of computation, we have reformu- 
lated the model of PHP in terms of classical Ising-spin- 
like variables, generalizing the corresponding formulation 
for the single-chain surfactants iQ. In this reformula- 
tion, a classical Ising-spin-like variable S is assigned to 
each lattice site; 5^ = 1 if the i-th lattice site is occu- 
pied by a water molecule. If the j-th site is occupied 
by a monomer belonging to a PHP then Sj = 1,-1 
depending on whether the monomer at the jth site is 
hydrophilic or, hydrophobic, respectively. The temper- 
ature T of the system is measured in the units of J/kB 
where J denotes the strength of the interaction between 
a spin and its six nearest-neighbours. This reformulation 
in terms of Ising-spin-like variables has been successfully 
used in studying a wide variety of phenomena exhibited 
by various types of surfactant molecules in aqueous me- 
dia pl- nsi a nd should not be confused with magnetic 
polymers||l9| . Besides, molecular dynamics simulation 
of similar molecular models pQ] have also been carried 
out to study the spontaneous formation of self-assemblies 
of surfactant molecules. 

Both the position of the center of mass and the confor- 
mation of the PHP is random in the initial state of the 
system. The allowed moves of the PHP are the same as 
those of the small surfactants in the Larson model (see ref 
]l2[ ) , namely, reptation, buckling and anti-buckling (also 
called pull) and kink movement |18|. Starting from the 
initial state, the system is allowed to evolve following the 



standard Metropolis algorithm: each of the attempts to 
move the PHP takes place certainly if A_E < and with 
a probability proportional to exp(— Ai<^/T) if AE > 0, 
where AE is the change in energy that would be caused 
by the proposed move of the PHP. 




In order to collect informations on the qualitative fea- 
tures of the conformations of the PHP we have directly 
looked at many snapshots of the PHP at various stages 
of MC updating of the state of the system. We have also 
computed several different quantities which provide im- 
portant quantitative informations on various aspects of 
the conformation of the PHP. 

A gross measure of the " size" of the PHP in water is 
given by its radius of gyration 
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where fj is the position vector of the j-th monomer and 
Rem is the position of the center of mass which is defined 

as Rem = (1/ia) J2j=i ^r 

Insight into the composition of the local neighbour- 
hood of an arbitrary hydrophilic monomer can be gained 
by computing the quantities Na, Nio and Ni^ which are 
the average numbers of its nearest-neighbour sites that 




FIG. 1. Typical snapshots of a PHP with La = 400 for 
(a) Lp = 4, (b) Lp = 40, (c) Lp = 50, and (d) Lp = 100. The 
hydrophilic and hydrophobic monomers are denoted by the 
black and grey circles, respectively. The boxes are unavoid- 
able artefacts of the graphics package used for plotting. 
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FIG. 2. The quantities Na, Nio,Niw, Noo and Now (see text 
for definitions) , which are represented collectively by the label 
A'', are plotted against / = Lp/La at a fixed temperature 
T = 2.0. The symbols +, x, *, open square and filled square 
correspond to Nii,Nio,Niw,Noo and A'^otu, respectively. 



are occupied by a hydrophilic monomer, a hydrophobic 
monomer and a water molecule, respectively. Similarly, 
the composition of the local neighbourhood of an arbi- 



trary hydrophobic monomer is reflected in the numbers 
Noi , Noo and Now , which are the average numbers of 
its nearest-neighbour sites that are occupied by a hy- 
drophiUc monomer, a hydrophobic monomer and a wa- 
ter molecule, respectively. Obviously, Nio = Noi as, 
throughout this paper, we consider PHP consisting of 
equal number of hydrophilic and hydrophobic segments 
of the same length Lp. 







FIG. 3. The non-zero elements of the contact maps of a 
PHP with La = 400 for (a) Lp = 4, (b) Lp = 40, (c) Lp = 50, 
and (d) Lp = 100 are denoted by dots. 



Suppose an index j {j — 1,2,..., La) labels the 
monomers sequentially along the primary structure of the 
PHP chain from one fixed end. The i, j-th element, dj, 
of the contact map C is defined to be non-zero if and 
only if in at least one of its equilibrium configurations 
the z-th and the j'-th monomers (irrespective of whether 
hydrophilic or hydrophobic) are not nearest-neighbours 
along the chain but occupy two nearest-neighbour lattice 
sites pi| . The contact map has been used to reconstruct 
the three-dimensional conformation of bio-polymcrs. 

For a given La, Lp and T , after equilibration, we have 
computed the above-mentioned quantities of our interest. 
Then we have repeated the calculations for several val- 
ues of La, Lp and T. All the data reported in this letter, 
however, have been generated for Lp = 400, correspond- 
ing to the longest PHP, for which we could sample, after 
equilibration, sufficiently large number of configurations 
required for averaging. 

For a fixed La — 400, typical snapshots of the PHP for 
a few different Lp are shown in the figs.la-d. The PHP 
is very stiff for Lp = A (fig. la). For intermediate values 
oi Lp, e.g., Lp = 40 (fig. lb) and Lp = 50 (fig.lc), it has a 
necklace- like conformation where "beads" of hydropho- 
bic monomers are connected by hydrophilic chains. Fi- 
nally, when Lp is of the same order as La, e.g., Lp = 100 
(fig. Id), the hydrophic monomers form a large collapsed 
globule surrounded by hydrophobic monomers. 

Each of the hydrophilic (hydrophobic) monomers has 
a tendency to have hydrophilic (hydrophobic) nearest- 
neighbours and avoid having hydrophobic (hydrophilic) 
nearest- neighbours. The snapshots shown in fig.l also 
indicate that a longer Lp enables the PHP to satisfy 
these tendencies. This can be shown more quantitatively 
(fig. 2) by plotting Nu, N^o,N^w,Noo and Now against 
/ = Lp/La at a fixed temperature T = 2.0. 
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FIG. 4. The average radius of gyration _R of a PHP is plot- 
ted against / — Lp/La at a fixed temperature T — 2.0. 

One striking feature of the PHP is that the shorter 
is the period the more stretched is the PHP, as shown 
by the snapshots in fig.l. This trend of variation is re- 
flected in the structure of the contact maps, shown in the 
flgs.3a-d, corresponding to the figs.la-d, respectively. In 
the contact map for / = Lp/La = 0.01 there are very 
few non-zero elements outside the diagonal backbone of 



the map. With increase of / more and more non-zero el- 
ements far from the diagonal backbone appear signalling 
folding or collapse of the PHP. This trend of variation of 
the " size" of the PHP can be seen quantitatively also in 
fig. 4 where we plot the radius of gyration R of the PHP 
as a function of /. 
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FIG. 5. The radius of gyration R is plotted against the 

temperature T at a fixed / = 0.01. 

Finally, keeping / fixed at a small value, say / — 0.01, 
corresponding to which the PHP is very stiff, if we raise 
T, the R of the PHP falls monotonically with increas- 
ing T (fig. 5), as expected, because of stronger thermal 
fluctuations. 

In summary, in this letter we have developed a Larson- 
type model of a periodic hetero-polymer. By carrying 
out MC simulations of these model PHP, each consist- 
ing of equal numbers of hydrophilic and hydrophobic 
monomers, we have investigated the effects of varying 
the period on its conformations in equilibrium. We have 
observed that, at a given temperature, the smaller is the 
ratio / = Lp/La, the stiffer is the PHP. We would like to 
emphasize that the stiffness of the PHP at a fixed temper- 
ature decreases with increasing Lp = La/{2n), where n is 
the number of segments of each type, in spite of the fact 
that, nLp, the total number of hydrophobic monomers 
remains fixed for the given La- This prediction, we be- 
lieve, can be tested directly in laboratory experiments. 

We thank E. Domany for enlightening discussions on 
contact maps and L. Santen for help with the graph- 
ics. This work is supported by the SFB341 Koln-Aachen- 
Jiilich and German-Israeli Foundation. 



[1] 

[2] 
[3] 
[4] 

[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 

[12] 

[13] 
[14] 
[15] 
[16] 
[17] 
[18] 

M. 

[20] 



B. Alberts, D. Bray, J. Lweis, M. Raff, K. Roberts 
and J.D. Watson, Molecular biology of the cell, (Gar- 
land Publishing, 1983); J. Darnell, H. Lodish and D. 
Baltimore, Molecular Cell Biology, (Scientific American 
Books, 1990) 

V.J. Pande, A. Yu Grossberg and T. Tanaka, Rev. Mod. 
Phys. (1999) 

B. Derrida, Phys. Rev. Lett. 45, 79 (1980); Phys. Rev. 
B 24, 2613 (1981). 

K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 
(1986); D. Chowdhury, Spin glasses and other frustrated 
systems, (Princeton University Press and World Scien- 
tific, 1986); K.H. Fischer and J. Hertz Spin Glasses 
(Cambridge University Press, 1991) 

E.I. Shakhnovich and A.M. Gutin, Nature 346, 773 
(1990) 

T. Garel, H. Orland and E. Pitard, in: Spin glasses and 
random fields, ed. A. P. Young (World Scientific, 1998) 
A.K. Chakraborty, E.I. Shakhnovich and V.S. Pande, J. 
Chem. Phys. 108, 1683 (1998) 

E. Orlandini and T. Garel, Eur. Phys. J. B 6, 101 (1998) 
M. Gerstein and M. Levitt, Sci. Amer. 279(5), 75 (1998) 
D. Stauffer, N. Jan and R.B. Pandey, Physica A, 198, 
401 (1993); D. Stauffer, N. Jan,Y. He, R.B. Pandey D.G. 
Marangoni and T. Smith-Palmer, J. Chem. Phys., 100, 
6934 (1994); N. Jan and D. Stauffer, J. Phys. II (France), 
4, 345 (1994). 

R.G. Larson, L.E. Scriven and H.T. Davis, J. Chem. 
Phys. 83, 2411 (1985); R.G. Larson, J. Chem. Phys. 89, 
1642 (1988); 91, 2479 (1989); J. Chem. Phys., 96, 7904 
(1992); Chem. Eng. Sci., 49, 2833 (1994). 
T.B. Liverpool, in: Annual Reviews of Computational 
Physics, vol. IV, edited by D. Stauffer, (World Scientific, 
Singapore 1996). 

F. Schmidt, in: Computational methods in colloid and 
interface science, ed. M. Borowko (Marcel Dekker, 1999) 
D. Stauffer and D. Woerman, J. de Physique II, 5, 1 
(1995). 

A.T. Bernardes, J. de Physique II, 6, 169 (1996); Lang- 

muir, 12, 5763 (1996). 

D. Chowdhury, J. de Physique II 5, 1469 (1995), Lang- 

muir, 12, 1098 (1996) 

D. Chowdhury, P.K. Maiti, S. Sabhapandit and P. 

Taneja, Phys. Rev. E 56, 667 (1997). 

P.K. Maiti and D. Chowdhury, Europhys. Lett. 41, 183 

(1998); J. Chem. Phys.109, 5126 (1998). 

Garel, H. Orland and E. Orlandini, cond 



T 



mat/9902147 



* On leave from the Physics Department, I.I.T., Kanpur 
208016, India. 



B. Smit, P.A.J. Hilbers and K. Essefink, Int. J. Mod. 
Phys. C 4, 393 (1993); S. Karaboni, K. Esselink, P.A.J. 
Hilbers, B. Smit, J. Karthauser, N.M. van Os and R. 
Zana, Science, 266, 254 (1994) 
[21] M. Vendruscolo, R. Najmanovich and E. Domany, Phys. 
Rev. Lett. 82, 656 (1999); L. Mirny and E. Domany, Pro- 
teins: Str., Func. and Genetics 26, 391 (1996); M. Vedr- 
uscolo, E. Kussell and E. Domany, Folding and Design 2, 
295 (1997) 



